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Abstract 

We study the exact dynamics underlying stimulated Raman adiabatic 
passage (STIRAP) for a particle in a multi-level anharmonic system (the 
infinite square- well) driven by two sequential laser pulses, each with con- 
stant carrier frequency. In phase space regions where the laser pulses 
create chaos, the particle can be transferred coherently into energy states 
different from those predicted by traditional STIRAP. It appears that 
a transition to chaos can provide a new tool to control the outcome of 
STIRAP. 



1 Introduction 

Laser radiation provides a means to control intra-molecular processes in a robust 
manner because of a conservation law that comes into play when monochromatic 
radiation interacts with nonlinear systems. The origin of this conservation law 
is the discrete time-translation invariance of laser driven systems. For radi- 
ation interacting with molecular systems, this conservation law gives rise to 
stable electron-photon (phonon-photon, roton-photon) structures described by 
conserved eigenstates (the Floquet states) of the driven system. Floquet states 
are exact eigenstates of periodically driven systems pQ, [2], 0- Their reality 
can be seen in recent atom-optic experiments 0], |S], where millions 

of sodium [S] or cesium atoms (4j, interacting with a time- modulated standing 
wave of light, underwent large coherent periodic oscillations in momentum in 
a multi-photon process. These coherent oscillations were subsequently found 
to be due to the interference of only two or three Floquet eigenstates whose 
phase space structure was determined by an underlying chaotic sea induced by 
the interaction between the atoms and the light 6 . In this paper, we wish to 
show that similar mechanisms are important when laser pulses interact with 
intra-molecular forces. We will focus on the exact dynamics underlying STI- 
RAP (stimulated Raman adiabatic passage) for laser pulses interacting with a 
simple model of intra-molecular dynamics. 

STIRAP has become an important tool for coherently controlling and chang- 
ing the vibration and electronic states of entire molecular populations with close 



1 



to 100% efficiency. STIRAP involves the application of short laser pulses with 
carefully chosen carrier frequencies to a molecular system for the purpose of 
exciting the molecules in a controlled manner. This technique causes a coherent 
change in the entire molecular population between molecular states. Traditional 
models of STIRAP generally view the molecular target as a simplified multilevel 
system. Indeed, traditional STIRAP focuses on three carefully chosen vibration 
and/or electronic levels. There are no studies, that we know of, that look at the 
effect on STIRAP of the actual nonlinear dynamics that occurs when the laser 
pulses interact with a molecular system. However, we know that laser pulses 
can induce chaos and this can strongly affect the response of the molecule to 
the laser field. In addition the internal dynamics of the molecule may itself be 
exhibiting the manifestations of chaos and simple pictures of the molecular level 
structure are likely not valid. Understanding this dynamics is very important 
for extending STIRAP to wider ranges of molecules 

In this paper we study the exact dynamics underlying STIRAP for a model 
system that contains the essential features of low energy vibration states, or 
rotational states, of a diatomic molecule driven by two laser pulses. We will find, 
for example, that in phase space regions where the laser pulses create chaos, the 
molecule can be transferred coherently into energy states very different from 
those predicted by traditional STIRAP. It appears that a transition to chaos 
may provide a new tool to control the outcome of these processes in molecular 
systems. 

STIRAP was first proposed by Hioe [7], [H| and later confirmed in an experi- 
ment involving population transfer between vibration-rotation states of sodium 
dimers [H], JH!- There are several variations to STIRAP which generally is de- 
scribed as a process involving three particular energy levels, Ei,E2, and E3, of a 
collection of atoms or molecules of interest. All atoms or molecules are initially 
in the lowest energy state Ei. Two laser pulses then impinge sequentially on 
the system in order to make a transition toward the target state E3. In the 
"ladder" version of STIRAP the target state £"3 is the highest state in energy 
and in the "lambda" version of the STIRAP the intermediate state E2 is the 
highest state in energy. The first laser pulse couples E2 and E3 which contain 
no population, and the second laser pulse then couples Ei and E2- As a result 
of these processes, the entire atomic or molecular population is transferred co- 
herently into the target state without populating the intermediate state after 
the laser pulses have passed [TT ] . P ] . P ] . [Hj. [Bj. [TH]. [T7 j . 

As mentioned above, conventional STIRAP analyzes an atomic or molecu- 
lar system in terms of three carefully chosen energy levels of the unperturbed 
system. However, complications arise in real applications since these systems 
are composed of multiple states due to the rotational, vibration, and electronic 
levels of the unperturbed system. Extension of STIRAP to multilevel systems 
has been discussed by a number of authors, generally in terms of simplified 
models using the rotating wave approximation and three or more laser pulses 
[IHllIlEniEIlESESIllESlEiEZIEHlESI- Alternatively, Raman chirped 
adiabatic passage schemes [23 I2H1 12S] > in which the frequency of the incident 
laser pulse is continuously changed, have also been proposed to allow the system 
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to climb through a sequence of molecular energy levels. 

In this paper we take a different approach to the problem. We consider 
the exact dynamics that takes place when the two laser pulses impinge on a 
multi-level system. Instead of isolating a certain number of levels of interest 
and analyzing the process in terms of the rotating wave approximation, we will 
let the full multilevel dynamics evolve and allow the system to tell us how many 
levels we must keep to accurately describe the atomic or molecular dynamics 
when a radiation field is present. Unlike the previous studies which utilize more 
than two pulses in multi-level systems, we only apply two pulses and allow 
as many levels as dictated by the dynamics participate in the process. One 
way to see how many levels must be kept is to look at the underlying classical 
phase space of the system of interest. The laser field, when interacting with 
nonlinear intra-molecular dynamics will induce nonlinear resonances and chaos 
in localized regions of the phase space. Those structures in the classical phase 
space which have a size greater than Planck's constant determine the structure 
of the Floquet eigenstates of the system, and thereby have a direct influence on 
the STIRAP process. 

The model we use to study the effect of chaos on STIRAP in a multi-level 
system is that of a particle in an infinite square- well potential |3l)j. |31| . The 
infinite square-well potential is an anharmonic potential of the form i^" in the 
limit n ^ oo. It can give some insight into the behavior of low lying vibration 
states, or rotation states, of molecular systems in the presence of sequential laser 
pulses. An approximate version of the square-well potential we consider here 
could also be constructed in an atom-optics experiment (Hj or in semi-conductor 
heterostructures. The dynamics of a particle in an infinite square-well potential 
is also interesting because the laser pulses can cause the low energy particle 
states to undergo a transition to chaos [221 1 123 • For the case of monochromatic 
laser fields, this transition to chaos is accompanied by a plateau of high harmonic 
radiation whose cut-off is determined by the width in energy of the chaotic sea 
induced by the laser field [5^ . 

In the sections below, we describe the behavior of a particle in an infinite 
square-well potential which is driven by two sequential laser pulses whose carrier 
frequencies are monochromatic and chosen to couple specific unperturbed energy 
levels of the particle in the square-well potential. We will find that for the case 
when the pulse amplitudes vary slowly in time so that the adiabatic theorem 
Pl] is satisfied, the dynamics can be described in terms "snapshots" of the 
underlying classical phase space at selected times as the laser pulses pass through 
the system. In Sect. 2, we describe the classical dynamics that results from the 
laser pulses. In Sect. 3, we discuss how we will describe the quantum dynamics 
for the driven system. In Sect. 4 we introduce Floquet theory. In Sects. 5, 
6 and 7 we show that a dynamics quite different from that of the traditional 
STIRAP ladder model can occur in multilevel systems due to the presence of 
chaos, even for fairly weak pulse amplitudes. Finally in Sect. 8 we make some 
concluding remarks. 
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2 Classical dynamics 



The model system for this study is a particle located in an infinite square well 
potential with spatial width, 2a. The potential energy has the form, V{x) = 
for |i| < a and V{x) — oo for > a, where x is the position of the parti- 
cle. The classical Hamiltonian which describes the dynamics of a particle of 
mass m moving in the potential well and driven by two sequential pulses of 
monochromatic radiation is given by 

H = h Uf{t)xcos{Cjft) + Us{t)xcos{Cjst), for |a;| < a, (1) 

where p is the momentum of the particle, t is the time, and ujf and Ug are the 
carrier frequencies of the first and second pulses, respectively. The amplitude 
of the first (second) pulse at time t is Uf{t) {Us{t)). If we rescale parameters 
using x = xa, p = ph/a, Uo = Uofi^ /2ma? , i = 2ma^t/h, H = HTt^ /2ma? , and 
Cjf^s — ujf\s^/2'ma'^, where h is Planck's constant, then the Hamiltonian is 

H ^ p'^ + Uf{t)xcos{ujft) + Us{t)xcos{uJst), for |a;| < 1 (2) 

and all parameters are dimensionless. The energy has been re-scaled in units of 
/2ma in order to make comparisons with the corresponding quantum system 
in subsequent sections. 

The amplitudes Uf{t) and Us{t) have Gaussian shape of the form, 

Uf{t) = t/oexp(-/3(t - tff) and Us{t) = t/oexp(-/3(t - Q^), (3) 

where tj < tg. We can control the duration of each pulse by adjusting the 
parameter (3 and we can control the amount of overlap of the two pulses by 
changing tf and ts- For simplicity, we assume that the maximum amplitude 
Uo and the width /3 of the two pulses are the same. A schematic picture of 
the variation in time of the amplitudes of the two pulses is displayed in Figure 
1. The first pulse is turned on and drives the particle in the square well and 
then, with an appropriate delay time, the second pulse is turned on. The whole 
pulse sequence ends at the total pulse duration time t = ttot- For the purpose 
of marking time intervals in our subsequent discussion, we choose times, ti = 
^ttot, tc — ^ttot, and ^2 — ^ttot- 

We will be interested in how the classical phase space behaves in the adiahatic 
limit where the pulses are turned on and off very slowly relative to certain 
intrinsic time scales in the system, such as the periods of the carrier frequencies. 
In this limit the amplitudes, Uf{t) and Us{t) remain essentially constant during 
time intervals where the cosine terms oscillate many times. We can get an idea 
of the structure of the phase space by plotting a Poincare surface of section 
of the phase space for fixed pulse amplitudes [H]. Thus we also consider the 
Hamiltonian 

H = p'^ + Uf{tfi^)xcos{ujft) + Us{tfix)xcos{LList), for |x| < 1, (4) 
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^ — -^cos((2z^- 1)6'- w^t), for 0<6'<7r. 



where the amplitudes are set to constants by choosing their value at some fixed 
time t — tjix- The Hamiltonian in Eq. ^is time-periodic and we can view the 
behavior of the phase space using Poincare surfaces of section [2] . The Poincare 
surfaces of section for time-periodic Hamiltonians are strobe plots of p and x, 

1. e. plots of p and x each time the Hamiltonian goes through one complete 
oscillation in time. 

We can perform a canonical transformation to action-angle variables (J, &) 
defined J = 2|p|/7r, and d — ±7r(a; -I- 1)/2. The Hamiltonian then has the form, 

H^^-^ E p^c„s((2.-l,.-.,0 

I/ — — OO ^ ' 

(5) 

An infinite number of nonlinear resonances are produced in the classical phase 
space by the external fields. The primary resonances are located ai J = = 
2a;/^s/((2j/ — l)7r^). As v increases, the energy at which higher order primary 
resonances are located decreases. 

In Figure 2, we show strobe plots of the classical phase space for the case with 
Uo = 3.0. The commensurability of these frequencies will allow us to use Floquet 
theory when we analyze the quantum system. We choose the pulse carrier fre- 
quencies to be a;/ = 37r^/4 and ojs — 57r^/4. For the frequencies we have chosen, 
the period of the Hamiltonian is To = 8/tt. For the five cases shown in Figs. 

2. a-2.e, we fix the amplitude of the pulses by setting (a) Uf^s{tfix)=Uf^s{ti), 
(b) Uf.Atf^.)=UfAtf), (c) UfAtf^.)=UfAtc), (d) UfAtf^.)=Uf.Ats), and (e) 
Uf.s{tfix)=Uf^s{t2), respectively. For each of these choices of amplitude Uf 
and Us, we show strobe plots of the classical phase space allowing the time 
dependence of the cosine waves to vary. The three largest primary resonances 
(zy — 1,2,3) induced by the first pulse are located at J = 2.5, J — 0.83 and 
J = 0.5, respectively. The three largest primary resonances {i' — 1, 2, 3) induced 
by the second pulse are located a,t J — 1.5, J — 0.5 and J = 0.3, respectively. In 
Fig. 2. a, with C//(ti) — 0.1667 and Us{ti) = 0.000003, the primary resonances 
induced by the first pulse are dominant. In Fig. 2.e, with Uf{t2) — 0.000003 
and Us{t2) — 0.1667, the primary resonances induced by the second pulse are 
dominant. In all cases, the first primary resonance [u = 1) is located at the 
highest energy and the higher order primary resonances are located at decreas- 
ing energy as v increases. As a result, this system will always have a chaotic 
region at low energy due to the overlap of higher order resonances. For energies 
above the region of influence of the u ~ I primary of the first pulse, the phase 
space is dominated by KAM (Kolmogorov-Arnold-Moser) tori. 

In Fig. 2.C, where tfi^ = ic = ■^ttot, the primary v — 1 resonances due 
to the two pulses have equal amplitude and are clearly visible at J — 2.5 and 
J = 1.5. For this case the pulse amplitudes are Uf{tc) — Us{tc) — 1.103. All 
the higher order primary resonances have been destroyed and a large chaotic 
sea has formed at low energy. 
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Fig. 2.b shows the classical phase space at time, tfi^ =tf when the first pulse 
reaches its maximum amplitude with Uf{tf) = 3.0 and Us{tf) = 0.055. The 
region of phase space about the primary v = 1 resonance due to the first pulse 
is very distorted by the resonance. Detailed calculation shows that there are 
small higher order (non-primary) resonance islands between J = 3 and J = 4. 
Fig. 2.d shows the classical phase space at time, tfix = tg when the second pulse 
has reached its maximum amplitude with U /{ts) = 0.055 and Us(ts) = 3.0. The 
primary v = 1 resonance due to the first pulse is very small and its region of 
influence does not extend very high in energy. 

It is interesting to compare the classical phase space for C/q = 3.0 with a case 
when the maximum pulse amplitude is [/o = 0.5. In Figure 3, we show the strobe 
plots of the classical phase space at the same times, tfix =t\, tf, tc, tg, t2 as 
in Fig. 2 but with Uo = 0.5. More island structures survive with this relatively 
weak value of pulse amplitude than in Fig. 2. The invariant surfaces between 
the two v = 1 primary resonances are distorted and higher order non-primary 
resonance islands can be seen even when the pulse amplitudes have reached 
their maximum values. 

3 The Quantum system 

The Schrodinger equation for the driven square-well system described in Sect. 
2 can be written (in dimensionless units) 



i — {x\^{t)) = -—^ + Uf{t)xcos{wft) + Us{t)xcos{wst) (a;|V(i)), (6) 



where the momentum operator is given hy p = —id/dx. In order to satisfy the 
boundary condition at x = ±1, the wavefunction should satisfy V'(a; = ±l,i) = 
{x = ±l|V'(t)) = for all times, t. 

For the case of a quantum particle in the infinite square well, when [//,s = 
(no driving force), the energy is conserved and we have a complete set of or- 
thonormal energy eigenstates which can be used as basis states to describe 
the dynamics of the driven system. For the unperturbed system, the en- 
ergy eigenvalues are En = n^7r^/4 and the orthonormal energy eigenstates are 
{x\En) = 4'n{x) = sin[n7r(a; — l)/2]. The dipole matrix elements in this basis 
are Xji 

= {En\x\En') where 



Note that integer values of J ( J = n) in the classical Hamiltonian correspond to 
the allowed quantized states of the quantum system. This simplifies comparison 
between the classical and quantum systems. 

We can expand the wavefunction, \il){t)), in the unperturbed energy basis so 
IV'(i)) = Tlin'^n{i)\En) ■ Then we can rewrite the Schrodinger equation in the 




[n + n'] (modulo 2) = 
[n + n'] (modulo 2) = 1. 



(7) 
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form 

= -iE„c„{t) + i[Uf{t)cos(ujft) + Usit)cos{uJst)] ^ it) (8) 

n' 

where Cn{t) = {En\tjj{t)) is the probabihty amplitude to find the system in the 
nth energy level at time t. We will generally assume that at time t = the 
system is in the state |V'(0)) = \Ei). We will then find the state \iIj{+oo)) after 
the two pulses have been turned on and off. 



4 Floquet States 

Once we fix the amplitudes, Uf{t = tfix) = Uf(tfix) and Us(t = tfix) = 
Us{tfix), the Hamiltonian becomes time periodic and the Schrodinger equation 
takes the form 

d / \ 

= [-^^ + Uf{tfix)xcos{ujft) + Us{tf,x)xco&{uJst)\ {x\%lj{t)). (9) 

For such systems, the energy is not conserved. However, if the carrier frequencies 
of the pulses are commensurate so w/ /ojs = /ng , where n / and are integers, 
then the Hamiltonian is invariant under a discrete time translation H(t) = 
H{t + To), where the period To of the Hamiltonian is 

To^J^ + ^). (10) 

For such systems, Floquet eigenstates, \4>ait)) (which have period To so \(j)a{t + 
To)) = \cf>a{t))) form a complete orthonormal basis which determines the dy- 
namics. Furthermore, the Floquet cigcnphases, fio., are conserved quantities 

We can obtain an eigenvalue equation relating fla and \(f)oi(t)). Consider 
the case when the system is in the ath Floquet eigenstate so that |'0(t)) = 
e~*^°'|(/)Q(t)). Then substitution into Eq. Q yields the eigenvalue equation 

H{t)-t-^yMt))=^Mt)), (11) 

where Hp{t)=H{t) — id/dt is the Floquet Hamiltonian. 

More generally, when the system is in the state \ip{0)) at time i = 0, the 
state of the system at time t can be written 

\m) = E Ao.e-^''"'\Mt)) - E {Mo)\mh-'''"'\Mt)), (12) 

a a 

The state of the system at time t = To takes on an especially simple form 

IV(ro)) =U;^(T„)|V(0)) =Ee-'''°^°l0a(O))(0.(O)|^(O)). (13) 
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where we have used the fact that \(f)a{To)) = |(/)q(0)). The Floquct evolution 
operator, \]p{To), can now be defined 

Vf(To) - ^e-^"°^°|</.„(O))(0„(O)|. (14) 

a 

We can compute matrix elements of the Floquet evolution operator in the basis 
of unperturbed energy eigenstates. Then the (n, n')th matrix element of the 
resulting Floquet matrix is given by 

Un,n'{To) = (S„|Uf(T„)|S„,) = ^e-'"°^<'(£;„|0,(O))(0„(O)|i;„,)- (15) 

a 

The ath eigenvalue of the Floquet matrix Un,n'{To) is exp(— ir^aTo) and the 
ath eigenvector in the unperturbed energy basis is given by a column matrix 
composed of matrix elements, (i?„|0Q(O)), where n — 1^ ...,oo. The eigenvalues 
fla can be obtained from exp(— iilo-To), but only modulus lOq- We refer to the 
eigenvalues Via obtained from the Floquet matrix as eigenphases. 

For the system we consider here, the Floquet matrix has a natural trunca- 
tion which is determined by the nonlinear dynamics of the system. Classically, 
the driven square-well system has a region of mixed phase space bounded at 
high energies by KAM tori. For the cases we will consider here, where the ini- 
tial state IV'(O)) is the unperturbed energy level |V'(0)) = the state \'ip{0)) 
can never penetrate very far into the high energy KAM region. This provides 
a natural truncation of the size of the Floquet matrix and we need to include 
only enough unperturbed basis states, \En), to cover adequately the region of 
mixed phase space. Each column of the Floquet matrix can be constructed by 
solving the time-dependent Schrodinger equation for one period. To, with the 
system initially in one of the unperturbed energy eigenstates. This integration 
is performed using each of the unperturbed energy eigenstates as an initial state 
until all the columns of the Floquet matrix have been computed. Floquet eigen- 
phases and eigenstates are obtained by numerically diagonalizing the Floquet 
matrix j^]. 

Husimi distributions allow us to visualize the distribution of probability of 
the Floquet eigenstates in the underlying classical phase space j35j. Physically 
they describe the location of the particle in the presence of the radiation field 
and provide important information about the actual dynamics taking place in 
the system. The Husimi distribution for a Floquet eigenstate \(j)a) is defined as 
H{xo,po) — \ {(l>a\xa,po)\'^ , where the state \xo,po) is a coherent state that can 
be represented in the position basis as 

/I \ f Y^"^ ( {x~XQf ipQ{x-XQ)\ 

(a;|a;o,po) = ( ^ I exp I + I. (16) 

The coherent state is a minimum uncertainty wave packet and has a parameter, 
cr, that determines the relative dispersion in both position and momentum space. 

In the subsequent sections, we will consider three different choices for carrier 
frequencies of the pairs of pulses which drive the system. For Case I, the first 
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pulse connects levels n = 2 and n = 3 and the second pulse connects levels n = 1 
and n = 2. This is the traditional model for the STIRAP ladder process 
However, as distinct from the usual discussion of STIRAP we will deal with the 
exact dynamics of the system. We will take account of the fact that we have a 
multilevel system that can undergo a transition to chaos. We will examine the 
effect of the full nonlinear dynamics on this system. For Case II, the first pulse 
connects levels n = 4 and n = 5 and the second pulse connects levels n = 1 
and n = 4. This is again a ladder process. For this case the underlying chaotic 
dynamics will have a surprising effect on the transition probabilities. Finally, 
for Case III, we consider a lambda process in which the first pulse connects levels 
n = 3 and n = 4 and the second pulse connects levels n = 1 and n = 4. In all 
cases, we consider the exact dynamics of the driven system. 

5 Case I: First pulse 2^3; Second pulse 1^2 

In this section, we examine the dynamics of the driven square- well system when 
two pulses are applied such that the first pulse connects levels n — 2 and n — 3 
and the second pulse then connects the levels n = 1 and n = 2. We first 
determine the behavior of Floquet eigenstates at fixed times t — tfix during 
which the pulses drive the system. The distribution of probability in the Floquet 
eigenstates is sensitive to structures in the classical phase space which are larger 
than Planck's constant. We then compare the prediction of Floquet theory to 
the actual behavior of the system in the non-adiabatic and adiabatic regimes. 

5.1 Behavior of Floquet Eigenstates 

The first pulse has carrier frequency cj/ = (£"3 — E2) = 57r^/4 and the second 
pulse has carrier frequency uJs — {E2 ~ Ei) = Stt^/A. These frequencies are 
commensurate since uf /us — 5/3. From Eq. Hl()|l . the period of the Hamiltonian 
is To = S/tt and the Floquet frequeney is Uo = Stt/Tq = 7r^/4. Thus, ujf = 5wo 
and UJs — 3ujo- We set Uo = 3.0. 

The dynamics of this system tells us that we only need to keep five unper- 
turbed energy eigenstates as a basis to form the Floquet matrix. This can be 
seen from Fig. 2 where we show the underlying classical phase space at selected 
values of tfix during the time that the pulses are on. For J > 5, the classical 
phase space is dominated by KAM tori with almost constant values of J and 
the unperturbed energy states are very weakly coupled by the dynamics for 
n > 5. Thus, to describe the quantum behavior of this system, it is sufficient 
to construct a 5x5 Floquet matrix with the five basis states \Ei), ...,\E5). We 
find that only four of the five eigenstates of the Floquet matrix are actively in- 
volved in the dynamics. Their eigenphases, fla, are plotted modulo ujo = 7r^/4 
in Figure 4. a. Two of these Floquet eigenphases are almost degenerate over 
the time interval that the pulses act and are not distinguishable on the scale 
shown in Fig. 4. a. For tfi^ — 0, the four Floquet eigenphases are approximately 
degenerate modulo ujo- 
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Wc can follow each Floquct cigcnstatc during the entire process by comput- 
ing the eigenstates for a sequence of values oitfix over the interval 0<tfix<ttot- 
For closely spaced values of t fix , Floquet eigenstates at different times belonging 
to different cigcnphascs will be orthogonal. This provides a means of following 
the evolution of each eigenstate as a function oft fix- As we will see, the Floquet 
eigenstates can change structure when avoided crossings occur between Floquet 
eigcnphascs. To keep track of the changes that occur in the Floquct eigenstates, 
we will give each eigenstate a unique alphabetical label determined by its dom- 
inant dependence on unperturbed energy states at time tfix = 0. We find that 
at tfix = the Floquct eigenstates have the following structure and we give 
them the following labels: 

A = |<^i) = D = \ct>i) = \E4,) a^nd E = \cj>^) = \E^) , 

BC+ = \<I>^) = ^{\E2) + \E^)), BC- = \<I>^) = ^{\E^)-\E^)). (17) 

The traditional STIRAP ladder process assumes that the molecule or atom 
in question can be approximated by a three level system and causes a coherent 
population shift of the atom from level 1 to level 3. We find that the traditional 
STIRAP ladder process occurs in our system for Uo < 0.1. For amplitudes 
Uo < 0.1, state D = 1^4) does not participate in the dynamics at all. The Flo- 
quet eigenvalue curve for ^^4, plotted as a function of tfix, crosses that for fii 
in two places but does not undergo any avoided crossings. The state D = \(l)4) 
remains predominantly dependent on \E4) during the entire process. For tradi- 
tional STIRAP, curve D in Fig. 4.b is replaced by state A and the only avoided 
crossing that occurs is the three-way avoided crossing at t fix = tc between states 
A, BC^ , and BC^ . The state A, which is predominantly composed of the state 
\Ei) before the multiple avoided crossing at t = tc, becomes predominately de- 
pendent on state |£^3) after the multiple avoided crossing, having interchanged 
its "1" character with the "3" character of states BC~^ and BC~ at the avoided 
crossing. Thus, at the end of the process state A would be composed predomi- 
nately of level n = 3 and the states BC+ and BC~ would be predominately of 
superpositions of levels n = I and n = 2. 

Once the amplitude Uo becomes greater than Uq = 0.1, something different 
happens due to the avoided crossing shown in Fig. 4.c. For Uo < 0.1 the Floquet 
eigcnphases for states A and D in Fig. 4.c cross at time tfix = Tj^^^ttot just 
before tfix = tc- For Uo > 0.1 the Floquet eigenphases for states A and D in 
Fig. 4.C avoid crossing at time tfix = t/. Before time tfix = t/ Floquet state A 
is predominately composed of level n = I and Floquct state D is predominately 
composed of level n = 4. After the avoided crossing at time t fix = tj the states 
have changed their character and Floquet state A is composed predominately 
of level n = 4 and Floquet state D is composed predominately of level n = I. 
Because of the avoided crossing at tfix = ti, the entire population gets shifted 
from level n = 1 to level n = 4 before the traditional STIRAP ladder process 
can take place. The traditional STIRAP ladder process now occurs among 
unpopulated states. It is interesting to note that the Floquet states A and D 
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cross at time tfix = r/7«13/23itot- A symmetry that was broken earlier appears 
to have been restored. These transitions are clearly seen in Fig. 5 where we 
show the level compositions of the four participating Floquet states, A, BC~^ , 
BC~ and Z? as a fimction of tfix- 

5.2 The Population Transfer 

Let us now determine the exact behavior of the system, when the pulses are 
applied, by solving the Schrodinger equation in Eq. |S1 We will assume that 
at time t = the system is in state |V'(0)) = l^'i)- As we will see, the actual 
dynamics of this system is determined by the length of time during which the 
pulses are allowed to act. The pulse duration time necessary to achieve adi- 
abatic behavior of the system is determined largely by the avoided crossings 
in the Floquet eigenphases. At isolated avoided crossings, involving only two 
eigenstates, the states involved interchange their character. 

Avoided crossings of Floquet eigenphases occur as the classical phase space 
becomes chaotic, and a symmetry has been broken in that local region of the 
phase space. The probability P^z that a transition occurs between the two 
Floquet eigenstates involved in an isolated avoided crossing can be computed 
from a formula obtained independently by Landau \6iS\ and Zener [57| . For our 
system, the Landau-Zener probability is given by 



where 5e is the eigenphase spacing at the avoided crossing and 7 is the rate of 
change of the Floquet eigenphases with respect to time t fix in the neighborhood 
of the avoided crossing. 

We have computed the Landau-Zener probability P^z for the isolated sharp 
avoided crossing at time tfix = ti shown in Fig. 4.c. The Landau-Zener prob- 
ability depends on ttot- The larger ttot^ the more "stretched out" the hori- 
zontal axis in Fig. 4.c will be relative to the vertical axis. We have obtained 
the following results by analyzing Fig. 4.c for different values of ttot- For 
ttot = 120, 6e = 0.005, and 7 ~ 0.001875 giving a Landau-Zener probability 
P^z = 0.979270. For ttot = 21000, 6e = 0.0063, and 7 = 0.00001376 giving a 
Landau-Zener probabihty Plz = 0.0108. For ttot = 270000, Se = 0.0030, and 
7 = 1.2x10^^ giving a Landau-Zener probability P^z^O- The first case is not 
in the adiabatic regime, but the second two cases are in the adiabatic regime 
because the probability of a transition is negligible. 

In Fig. 6, we show the probability — \{E„\ip{t))\'^ (for the four levels 

n = 1, 2, 3, 4) to find the system in the nth unperturbed level at time t for 
the three cases; ttot — 120, ttot = 21000 and ttot = 270000. These results are 
obtained by directly solving the Schrodinger equation, Eq. ||SJ). In all cases 
we start the system in the initial state, ItpiO)) = l^-i) with maximum pulse 
strength, Uo = 3.0. In Fig. 6. a, where there is a large Landau-Zener probability 
for the system to jump from Floquet state A to Floquet state D, the system 




(18) 
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comes out of the sharp avoided crossing at tfix = tj still predominately in the 
level \Ei), and the traditional STIRAP ladder process can then occur at i = <c- 
As the pulses are turned on and off, the system transitions from level \Ei) to 
level \Es). In Figs. 6.b and 6.c, the Landau-Zener probability is essentially zero 
and no transition occurs at the sharp avoided crossing at tfi^ = ti. The system 
comes out of the sharp avoided crossing in level I-E4). As the laser pulses are 
turned on and off the system transitions from the initial state |V'(0)) = \Ei) 
to the final state |V'(+oo)) = Note that both Figs. 6.b and 6.c follow 

almost exactly the behavior of the Floquet state A shown in Fig. 5.c. This is 
an indication that we are in the adiabatic regime in Figs. 6.b and 6.c. 

The very large oscillations in the probability in Figs. 6.b and 6.c have been 
explained by Berry |2H| in terms of a sequence of " super-adiabatic bases" . He 
shows that the decrease in the amplitude of these oscillations as we increase ttot 
is a sign that we are moving further into the adiabatic regime. The frequencies 
of the oscillations in Figs. 6.b and 6.c appear to be determined by the difference 
in Floquet eigenphases of the two Floquet states involved in the sharp avoided 
crossing. For example, at t fix = tf the period of the oscillation is rosc~400. The 
difference in the Floquet eigenphases is |AJ7| = — ri4|ft!0.016. Thus, Tosc = 
27r/|Arj| = 393. Similarly, at tfi^ = (r/-i/)/2, Tosc~GOO. The difference in the 
Floquet eigenphases is \An\ = {n^ - O4|«0.011. Thus, Tosc = 27r/|AO| = 571. 
The observed oscillation periods are the same for both Fig. 6.b and Fig. 6.c. 

6 Case II: First pulse 4^5; Second pulse 1^4 

We now turn on pulses with higher carrier frequencies in order to examine more 
closely the relation between the quantum transitions and their relation to the 
underlying classical dynamics. We first apply a pulse whose carrier frequency is 
ujf = {E5 — E4) ~ 9tt^/A. We then apply a second pulse whose carrier frequency 
is LUs = {E4 — El) — 157r^/4. The two frequencies are commensurate since 
ojf/ujs = 3/5. From Eq. (|10|l . the period of the Hamiltonian and the Floquet 
frequency are again To = S/tt and ujo = 7r^/4, respectively. Thus, w/ = 9uJo and 
LOs = 15a;o. We will consider the case when the maximum amplitude of both 
pulses is Uo = 13.0. For these frequencies and amplitudes, we find that we can 
induce a transition of the entire population from level \Ei) to level \Eiq). Below 
we describe how this happens. 

Classical phase space plots for times tfi^ — tf, tfi^ — tc, tfix — tjv = 
^/^ttot, and tfix — ts are shown in Figs. 7.a-7.d, respectively. The first primary 
resonance due to the first pulse is located at J = 4.5 and the first primary 
resonance due to the second pulse is located at J = 7.5. The frequency of the 
second pulse is chosen to connect levels n = 1 and n = 4. However it also 
connects the levels n = 7 and n = 8 since Eg — E-j — Idujo- This is why the first 
primary resonance due to the second pulse lies at J = 7.5. The states below 
J = 4.0 are immersed in the chaotic sea formed by the higher order primaries 
induced by the two pulses during most of the time that the one or the other of the 
pulses have a significant strength. Higher-order non-primary resonance islands 
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can be seen above J — 6.0 during most of the pulse sequence. In Fig. 7.b there 
is a chaotic sea which has formed throughout the region J = to J = 9. Also, 
during the last half of the pulse sequence, the KAM tori near J = 10 are highly 
distorted due to the formation of the primary resonance at J = 7.5. Thus, from 
the classical phase space we see that it requires approximately twelve square- well 
energy eigenstates to accurately describe the dynamics of this system. 

The Floquet matrix that we use to describe the quantum dynamics is a 
12x12 matrix. However, we find that only ten Floquet eigenstates are directly 
involved in the dynamics. To keep track of these ten Floquet eigenstates, we 
will give each state a unique alphabetical label determined by their dominant 
dependence on unperturbed energy states at time tfix = 0. We find that at 
tfix — the Floquet eigenstates have the following structure and we give them 
the following labels: 



A=|0l) = |Sl>, B=|02> = |S2>, |03> - 1^3), i^=|06>-|S6> 

DE+ = |04> = ^(|i?4> + \E,)), DE- = |05> - -^m) - |i?5», 

G= |<^7) = I^^t), = 108) = 1^8), 109) = 1^9), J= l^io) = I^IO). 



(19) 



The Floquet eigenphases corresponding to these ten Floquet eigenstates are 
plotted modulo lOq in Fig. 8. A number of avoided crossings occur between 
the eigenphases during the time the pulses act on the system. There are four 
avoided crossings that largely determine the dynamics. There is a multiple wide 
avoided crossing at tfi^ = tc = ^ttot which involves the seven states, B, C, 
DE^, F, H, and /. There is a sharp avoided crossing at ij^x — 'Tiii^'ottot that 
involves the states B and G. There is a three-state wide avoided crossing at 
tfix — Tiv^^ttot that involves states A, B and H. There is a sharp avoided 
crossing at tfix = Tv^^ttot which involves the states A and J. 

In Fig. 9. a we plot the eigenphases of the seven Floquet states B, C, DE^, 
F, H, and / involved in the multiple avoided crossing at tfix = tc = \ttot- 
These states have support in the unperturbed square- well levels -Es, -B4, -E5, Eq, 
i?7, i?8, and Eq. In Fig. 9.b we show a magnification of the very sharp avoided 
crossing between states DE^ and H at tf ix — tc — \ttot- The effect of these 
avoided crossings can be seen in the dependence of the Floquet eigenstates on 
the square- well states |-E„). Plots of Floquet eigenstates B, H, and G are shown 
in Fig. 10 and plots of Floquet eigenstates DE^, G and F are shown in Fig. 11. 
There is a complicated interchange of levels occurring. As shown in 33 , at such 
multiple avoided crossings, the Floquet eigenstates emerge with very different 
probability distributions than the entering states. Thus, multiple avoided cross- 
ings behave differently from isolated pairs of avoided crossings where the states 
simply interchange character. Multiple avoided crossings provide a mechanism 
for the spread of the manifestations of chaos in quantum systems |33j . 

It is useful to note that there is an isolated avoided crossing at tfix = tuj 
that causes states B and G to switch from i?«|i?2) and Cwji^y), as they enter 
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the avoided crossing at tju, to B^lEr) and Gk\E2) as they leave. Thus state 
BKi\Ej) as it enters the multiple avoided crossing at tfi^ — tc — ^Uot- 

Let us now consider the transition that causes the population of the square- 
well to undergo a coherent transition from level \Ei) to level jE'io)- This can 
occur if the system evolves adiabatically and follows the behavior of Floquet 
state A. The level dependence of Floquet state A is shown in Fig. 12. a. It 
starts out in level \Ei) and then partially switches to level \E^) at tfix = tjv 
due to a three state avoided crossing between states A, B, and H, and finally at 
tfix = TV it switches completely to level \Ew) due to a sharp avoided crossing 
between states A and J. The avoided crossings at times t fix = tiv and t fix = ry 
that cause these transitions are shown in Figs. 13. a and 13. b, respectively. 

In Fig. 14 we show a sequence of Husimi plots of the states B, H, A, and 
J as they go through the avoided crossings at tfix = tjv and tfix — Ty- The 
Husimi plots show the location of the quantum particle in phase space when 
the system is in a given eigenstate. The columns from left to right, show the 
Floquet states B, H, A and J. From top to bottom, the Husimi plot of each 
state is shown at times (a) tfix = tf, (b) tfix = tc, (c) tfix = tiv, (d) tfix = Ty 
and (e) t fix — ^2 ■ The sequence of events causing Floquet state A to undergo a 
transition from level \Ei) to level \Eio) can be seen clearly in these plots. 

Finally in Figs. 15. a and 15.b, we plot the probability P„(t) = \{En\ip{t))\'^ 
(for the ten levels n = 1, 10) to find the system in the nth unperturbed level 
at time t (Fig. 15. c will be discussed in Sect. 7). The system is prepared in 
initial state |V'(0)) = l^^i) with maximum pulse strength, Uo = 13.0. We obtain 
these plots by directly solving the Schrodinger equation in Eq. JHIl. In Fig. 15. a 
we show the results for the nonadiabatic case. We choose ttot ~ 600 which is 
a rapid evolution of the pulses. After the pulses have passed, the final state of 
the system is \ip{+oo)) = lE^). The adiabatic case is shown in Fig. 15. b. We 
now choose ttot = 6000. The probabilities now closely follow the behavior of 
the single Floquet eigenstate A as can be seen in Fig. 12. a. After the pulses 
have passed, the final state of the system is |V'(+oo)) = \Eiq). This transition 
dynamics is determined by the structure of the avoided crossings in the Floquet 
eigenphase curves and by the phase space structure of the Floquet eigenstates 
at the avoided crossing. 

The phase space structure of the Floquet eigenstates (see the Husimi plots) 
is determined by structures in the underlying classical phase space which are 
greater than Planck's constant. The transition from |'0(O)) = l-^i) to \tp{oo)) = 
\Eio) occurs because at the time of the transition the classical phase space is 
connected by a chaotic sea from J = 1 to J = 9 and a strongly distorted KAM 
region at J = 10 which allows the quantum state to tunnel into the phase space 
region around J = 10. This can be seen from the Husimi plot for the state A 
at time tfix = |§^tot which is shown in Fig. 16. The unperturbed levels, n — 1, 
n — 5 and n = 10 are connected via the chaotic sea around the two v — 1 
primary resonances which exist at the time tfix = ^f^tot- 
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7 Case III: First pulse 3—^4; Second pulse 1—^4 

The final case we consider is a lambda process in which the system is driven by 
pulses with carrier frequencies w/ = {E4 — E^) = 77r^/4 and Ug = {E4 — Ei) = 
157r^/4. These frequencies arc commensurate and the period of the Hamiltonian 
and the Floquet frequency are again Tg = S/tt and ujo = 7r^/4, respectively. 
To stay as close to Case II parameters as possible, wc chose pulse amplitudes 
Uo = 13.0. Although this case looks similar to Case II, we will find quite different 
results. 

In the adiabatic limit, the transition \Ei)-^\E^) occurs, as shown in Fig. 
15. c, which is very different from Case II. We believe that the difference between 
these two processes is due to the difference in the phase space; distribution of the 
resonances and chaos induced in the system by the laser pulses. Strobe plots 
of the classical phase space for the lambda case for times tjix = tj, tfix = tc, 
tfix = Tjv = 0.62<;tot, and tfix = tg arc shown in Figs. 17.a-17.d, respectively. 
The first primary resonance due to the first pulse is located at J = 3.5 (rather 
than J = 4.5 as is in Case II) and the first primary resonance due to the second 
pulse is located at J = 7.5 (the same as Case II). For Case II, a pathway 
is opened by a chaotic sea that allows a Floquet state to tunnel across the 
entire energy region (see Fig. 16) from \Ei) to \Eiq) at the time of the avoided 
crossing. In the strobe plots of the classical phase space in Fig. 17 for Case III, 
this pathway is blocked by KAM tori which do not allow the quantum system 
to tunnel out of the low energy region. We also find fewer avoided crossings in 
the plot of the Floquet eigenphases for this case. 

8 Conclusion 

We have found that sequential laser pulses, when applied to a particle in an 
anharmonic multilevel system, can induce nonlinear resonances and a transition 
to chaos in the dynamics of the particle. The extent of the region influenced by 
resonances and chaos determines the number of unperturbed energy eigenstates 
that must be kept to form the basis used to construct the Floquet matrix and 
to determine the quantum dynamics. 

When the pulses are applied in a manner which adiabatically changes the 
dynamics of the particle, this transition to chaos can be used to control the co- 
herent transfer of the particle across the chaotic sea from a low lying energy state 
to a highly excited energy state. Floquet theory provides an accurate means 
of describing the dynamical behavior of these driven systems in the adiabatic 
limit. The type of population transfer that is allowed depends on the nature of 
the Floquet eigenphase avoided crossings created by the underlying transition to 
chaos and the ability of the corresponding Floquet eigenstates to tunnel across 
vast regions of the phase space because of the induced chaos. In the Cases I and 
II, avoided crossings occurred between Floquet eigenstates which had spread 
throughout the available phase space, and a coherent flow of probability from 
the low energy side to the high energy side of the chaotic sea occurred. In Case 
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Ill, the path to the high energy states around the resonance at J = 7.5 appears 
to be blocked by KAM tori and the transition to the higher energy states does 
not occur. 

In molecular systems, the success of STIRAP will depend on the state of 
the dynamics of the unperturbed molecular system, and on the changes in that 
dynamics induced by the laser pulses. Most molecular systems have regimes of 
internal chaos and the interplay of these regimes with the dynamics induced by 
the laser pulses is critical to understanding STIRAP in these system. 
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Figure 1: Schematic diagram for the two pulses. The first pulse connecting 
levels \E2) and ji^a) is shown as a solid line. The second pulse connecting levels 
\Ex) and I-E2) is shown as a dotted line. They have maximum strength Uo at 
times t = tf and t = tg, respectively. The whole pulse sequence takes a time 
t = ttot to complete. In the figure, ti = l/20ttot, tc = ^/^Uot and t2 = 19/20ttot- 

Figure 2: Strobe plots of the action- angle variables (J, 0) for the infinite square 
well system with pulse amplitudes Uo = 3.0 and frequencies uif = 5uio and 
LOs = 5lOo- Strobe plots arc shown at times (a) tji-j. = ti, (b) tfi^ — ^fi (c) 
tfix = tc, (d) tfix = tg and (e) tfix = t^- For each plot O<0<7r. The three 
largest primary resonances = 1, 2, 3 due to the first pulse are located at 
J = 2.5, 0.83 and 0.5. The three largest primary resonances^ = 1, 2, 3 due to 
the second pulse are located at J = 1.5, 0.5 and 0.3. 

Figure 3: Strobe plots of the action-angle variables (J, 9) for the infinite square 
well system with pulse amplitudes Uq = 0.5 and frequencies ujf = 5wo and 
u>s = 'dojo, respectively. Strobe plots are shown at times (a) tfix = ti, (b) 
tfix = tf, (c) tfix = tc, (d) tfix = ts and (e) tfix = t2- For each plot 0<6'<7r. 

Figure 4: Floquet eigenphases, for the system with maximum pulse strength 
Uo = 3.0 and frequencies w/ = 5wo and LOg = 'SlOo, arc plotted over the entire 
interval 0<tfix<ttot- (a) Floquet eigenphases for four Floquet states A, BC^ , 
BC~ and D plotted mod ujo = 7r^/4. (b) Floquet eigenphase curves for the 
Floquet states D, i?C+and BC^ . The three-level wide avoided crossing at 
tfix = tc is clear, (c) Floquet eigenphase curves for the Floquet states A and 
D. The sharp avoided crossing at tfix = ti and the crossing at tfix = th are 
clearly seen. 

Figure 5: Probability distribution |(i?„|(/)Q,)p of the unperturbed energy levels 
\En) which compose each of the Floquet eigcnstatcs (a) BC^ (b) BC^ , (c) 
A and (d) D, plotted over the entire interval 0<tfix<ttot for pulse amplitude 
Uo = 3.0 and frequencies w/ = 5uJo and ujg = Suo- The probability curve for 
level \En) is labeled with level quantum number n. 

Figure 6: The probability Pn{t) = |(£^„|V'(i))P to find the system in the un- 
perturbed level \En) for the system prepared in initial state \ip{0)) = l-E'i) with 
maximum pulse strength, Uo = 3.0 and frequencies w/ = 5u>o and Ug = 3uio- 
The total pulse duration times are (a) ttot = 120, (b) ttot = 21000, and (c) 
ttot = 270000. The numbers attached to each curve show the components of the 
transition probability in terms of the unperturbed energy eigenstate basis. Case 
(a) is not in the adiabatic regime. Cases (b) and (c) are within the adiabatic 
regime and basically reproduce the structure of the single Floquet eigenstate A 
in Fig. 5.C. 
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Figure 7: Strobe plots of the action-angle variables (J, 6) for the infinite square 
well system with pulse amplitudes Uo — 13.0 and frequencies LOf = and 
ojs = 15ciJo, respectively. Strobe plots are shown at times (a) tfix = i/, (b) 
tfix = ic^ (c) tfix = Tiv = i/bttot and (d) tfix = tg. The first primary resonance 
from the first pulse is located at J = 4.5 and the first primary resonance from 
the second pulse is located at J = 7.5. 

Figure 8: The ten Floquet eigenphases, plotted modulo Wo, which determine the 
dynamics for pulse strength Uo — 13.0 and frequencies ujf = and ojs — 15cJo- 
The curves are identified following the classification scheme in Eq. [T^ 

Figure 9: (a) The seven Floquet eigenphases, plotted modulo loo, which are 
involved in the multiple avoided crossing at tfix = tc = \ttot for pulse strength 
Uo = 13.0 and frequencies Wf = 9aJo and ujg — IbuJo- The curves are identified 
following the classification scheme in Eq. E| (b) A magnification of the very 
sharp avoided crossing between Floquet states DE^ and H at tfix — tc — ^ttot- 

Figure 10: Probability distribution |(i?„|(/)a)p of the unperturbed energy levels 
which compose each of the Floquet eigenstates (a) B, (b) H and (c) G, plotted 
over the entire interval 0<tfix<ttot for pulse strength Uo = 13.0 and frequencies 
cof — 9uJo and cOs = 15wo- The probability curve for level |i?„) is labeled with 
level quantum number n. 

Figure 11: Probability distribution Kii'nl^Q)^ of the unperturbed energy levels 
which compose each of the Floquet eigenstates (a) DE~ , (b) DE~^ , (c) C and 
(d) F plotted over the entire interval 0<t fix<ttot for pulse strength Uo — 13.0 
and frequencies uif — 9ujo and ajg — IButo- The probability curve for level \En) 
is labeled with level quantum number n. 

Figure 12: Probability distribution |(i?„|0Q)p of the unperturbed energy levels 
which compose each of the Floquet eigenstates (a) A and (b) J, plotted over 
the entire interval 0<tfix<ttot for pulse strength Uo = 13.0 and frequencies 
ujf = 9uJo and ajg = i^ujo- The probability curve for level \En) is labeled with 
level quantum number n. 

Figure 13: (a) Magnification of the three Floquet eigenphases for eigenstates 
A, B and H which are involved in the three-state avoided crossing at tfix = 
Tiv = i/bttot for pulse strength Uo — 13.0 and frequencies tof — 9lOo and 
LOg = XbiJo- The curves are identified following the classification scheme in Eq. 
1191 (b) Magnification of the sharp avoided crossing at tfix = Ty = 2/ittot for 
pulse strength Uo — 13.0 and frequencies ojf = 9ujo and LOg = 15wo. Curves of 
eigenphases for Floquet eigenstates A and J are shown. 



20 



Figure 14: Husimi plots for the Floquet states (a) B, (b) H, (c) A and (d) J 
plotted in columns from left to right, respectively, for maximum pulse strength 
Uo ~ 13.0 and frequencies LOf = 9ujo and LUg = 15lUo- Each Floquet state 
(column) is shown (from top to bottom) at times (a) tfix = tf, (b) tfix = tc, 
(c) tfix = Tiv = 3/5ttot, (d) tfix = TV = 2/3ttot, (e) tfix = t2- The effect of 
the avoided crossings that enable Floquet state A to undergo a transition from 
level \Ei) to level |£'io) can be seen clearly. 



Figure 15: The probability |(_E„|V'(t))P to find the system in the unperturbed 
level \En) for the system prepared in initial state |V'(0)) = \Ei) with maximum 
pulse strength, Uo — 13.0 for the ladder processes in (a) and (b) and for the 
lambda process (c). For the ladder process, the frequencies of both pulses are 
tOf = 9uJo and lJs ~ IbWo and (a) ttot = 600 and (b) ttot = 6000. For the lambda 
process, the frequencies of both pulses are Wf = IlOq and uog = ISWo and (c) 
ttot = 6000. The numbers attached to each curve show the components of the 
transition probability in terms of the unperturbed energy levels. 



Figure 16: Husimi plot for the Floquet state A for maximum pulse strength 
Uo = 13.0 and frequencies ojf = 9uJo and ajg = 15ujo at time tfix = 13/20ttoi5 
just before the avoided crossing at tfix = Ty- This state covers the entire 
chaotic region of the underlying classical phase space from J = 1 to J = 9 and 
the highly distorted mixed region at J = 10. 



Figure 17: Strobe plots of the action-angle variables (J, 6) for the infinite square 
well system with pulse pulse amplitudes Uo = 13.0 and frequencies Wf = TuJo and 
ujs — 15ujo, respectively. Strobe plots are shown at times (a) tfix =tf, (b) tfix = 
tc, (c) tfix = Tiv = 0.62ttot and (d) tfix = ts- The first primary resonance from 
the first pulse is located at J = 3.5 and the first primary resonance from the 
second pulse is located at J = 7.5. 
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